BCS-BEC crossover in 3D Fermi gases with spherical spin-orbit coupling 
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We present a systematic theoretical study of the BCS-BEC crossover problem in 3D atomic Fermi gases at 
zero temperature with a spherical spin-orbit coupling which can be realized by a synthetic non-Abelian gauge 
field coupled to fermions. Our investigations are based on the path integral formalism which is a powerful 
theoretical scheme for the study of the properties of bound state, the superfluid ground state, and the collective 
excitations in the BCS-BEC crossover. At large spin-orbit coupling, the system enters the BEC state of a novel 
type of bound state (referred to as rashbon) which possesses a non-trivial effective mass. Analytical results and 
interesting universal behaviors for various physical quantities at large spin-orbit coupling are obtained. Our 
theoretical predictions can be tested in future experiments of cold Fermi gases with 3D spherical spin-orbit 
coupling. 
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I. INTRODUCTION 

It has been widely accepted for a long time that, by tuning 
the attractive strength in a Fermi gas, one can realize a smooth 
crossover from the Bardeen-Cooper-Schrieffer (BCS) su- 
perfluidity at weak attraction to Bose-Einstein condensation 
(BEC) of difermion molecules at strong attraction lll-j^. For 
a dilute Fermi gas in three dimensions where the effective 
range ro of the short range interaction is much smaller than 
the inter-particle distance, the system can be characterized by 
a dimensionless parameter l/(^Faj), where a, is the s-wave 
scattering length of the short range interaction and fcp is the 
Fermi momentum in the absence of interaction. The BCS- 
BEC crossover occurs when the parameter l/(^Ffl.s) is tuned 
from negative to positive values, and the BCS and BEC limits 
correspond to the cases llik^as) -oo and l/ikfOs) -t-oo, 
respectively. The BCS-BEC crossover has also become an 
interesting issue for the studies of dense nuclear a nd q uark 
matter which may exists in the core of compact stars |[lC)l - [l2ll . 

This BCS-BEC crossover phenomenon has been success- 
fully demonstrated in ultracold fermionic atoms, where the s- 
wave scattering length and hence the parameter 1 / (kpOs) was 
tuned by means of the Feshbach resonance lfl3l - [l5ll . At the 
resonant point or the so-called unitary point where Us — > oo, 
the only length scale of the system is the inter-particle dis- 
tance (~ fcp'). Therefore, the properties of the system at the 
unitary point l/{kpas) - become universal, i.e., indepen- 
dent of the details of the interactions. All physical quanti- 
ties, scaled by their counterparts for the non-interacting Fermi 
gases, become universal constants. Determining these univer- 
sal constants has been one of the most intriguing topics in the 
research of the cold Fermi gases lfl6l[l7ll . 

While the BCS-BEC crossover triggered by tuning the 
attraction strength between fermions from weak to strong 
(l/ikpOs) from -oo to -i-oo) has been comprehensively stud- 
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ied both theoretically and experimentally, it is always inter- 
esting to look for other mechanisms to realize the BCS-BEC 
crossover. Recent experimental breakthrough in generating 
synthetic non-Abelian gauge field has opened up the oppor- 
tunity to study the spin-orbit coupling (SOC) effect in cold 
atomic gases lITsll . For fermionic atoms fl^, it may provide 
an alternative way to realize the BCS-BEC crossover 112011 . 

Such SOC for fermions can be realized by generating a 
synthetic SU(2) gauge field A'' of the type = Ai^ U^. 
From the minimum coupling scheme, the resulting Hamilto- 
nian for fermions reads - p-/(2m) - t ■ where - 
(A^p_x, Ayfy, A,p^). The gauge field strengths Ai [i - x,y,z) 
characterize the spin-orbital coupling constants. The problem 
of the difermion bound state in 3D case in the presence of 
SOC has been studied in ll2lll . Three special cases were con- 
sidered: (!) Ajf - Ay - and A^ - A (called extreme prolate 
(EP)); (2) A;c = Ay - A and A, - (called extreme oblate 
(EO)); (3) A_, = Ay ^ A, ^ A (called spherical (S)). The EO- 
type SOC is physically equivalent to the Rashba SOC which 
is interesting for condensed matter physics. For EO- and S- 
type SOCs, it was shown that the difermion bound state exists 
even for a, < where the bound state does not exist in the 
absence of SOC. With increased SOC, the bounding energy 
is generally enhanced[21]. The bound state also possesses a 
non-trivial effective mass which is generally larger than twice 
of the fermion mass m ll22l - l24ll . Such a novel bound state 
caused by the SOC is now referred to as rashbons in the liter- 
atures 122]. For the 2D case, the bound state exists for arbi- 
trarily small attraction. It was shown in ll25ll that the EO-type 
SOC or the Rashba SOC enhances the binding energy and the 
bound state also has a non-trivial effective mass. This is anal- 
ogous to the catalysis of the dynamical mass generation by an 
external non-Abelian gauge field in quantum field theory ll26ll . 

Because of the presence of novel bound state with SOC, it 
has been proposed that a dilute Fermi gas with EO- and S- 
type SOCs can undergo a smooth crossover from the BCS 
superfluid state to the Bose-Einstein condensation of rash- 
bons (RBEC) even for negative values of l/(A:Ffl.5) if the 
SOC constant A is tuned from small to large values ll20ll . 
Due to the presence of SOC constant A, the 3D BCS-BEC 
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crossover problem depends on two dimensionless parameters: 
l/ikfOs) and A/kp (we set m = 1 in this paper). The BCS- 
BEC crossover problem and anisotropic superfluidity in 3D 
Ferr ni g ases with EO-type SOC has been extensively stud- 
ied I23I I27I1 . It was shown that the system enters the RBEC 
regime at A/kp ~ 1 for EO-type SOC for negative values of 
l/ikpOs). The BCS-BEC crossover in 2D Fermi gases with 
EO-type SOC was also studied i25i,28i1. Similar conclusions 
were found for the 2D case. 

In this paper, we present a systematic theoretical study of 
the BCS-BEC crossover in 3D Fermi gases at zero tempera- 
ture with S-type SOC. Especially, we will study the properties 
of the collective modes along the BCS-BEC crossover and the 
effective interaction among the rashbons in the RBEC regime. 
As far as we know, in the presence of SOC, these two interest- 
ing issues has not yet been studied (See the note added). For 
S-type SOC, the superfluid ground state is isotropic, which 
brings much convenience to the computations, and enables us 
to obtain various analytical results and universal behaviors at 
large SOC. 

This paper is organized as follows. In Section II, we set 
up the functional path integral formalism for the BCS-BEC 
crossover problem with a spherical SOC. Then we first deter- 
mine the binding energy and the effective mass of the rashbon 
at vanishing density and temperature (the vacuum in the pres- 
ence of SOC) in Section III. The ground state properties, such 
as the solution of the gap and number equations, fermion mo- 
mentum distribution, the condensate fraction, and the super- 
fluid density are discussed in Section IV. We derive the Gross- 
Pitaevskii free energy for the weakly interacting rashbon con- 
densate at large SOC and determine the rashbon-rashbon scat- 
tering length in Section V. The properties of the collective ex- 
citations, such as the gapless Goldstone mode and the mas- 
sive Anderson-Higgs mode, are investigated in Section VI. We 
summarize in Section VII. 



II. MODEL AND EFFECTIVE POTENTIAL 

The Hamiltonian describing spin- 1/2 fermions moving in 
three spatial dimensions in a SU(2) non-Abelian gauge field 
A^ is given by 



HcF = I dhif,\r) 



,(P -A%) 

2m 



«A(r) 



(1) 



where t//{r) = i/'|(r)]^ represents the two-component 

fermion fields and p = -/SV is the momentum operator. Here 
are Pauli spin operators (// - x,y,z)- In the following we 
use the natural units ti - k^ — m — \. Considering a spherical 
uniform gauge field, = -A5^^ {i = x, y, z), the Hamiltonian 
can be reduced to 

HoF = j d\ilf\r) + ■ pj iA(r), (2) 

which can be regarded as a generalized Rashba Hamilto- 
nian with an isotropic SOC. Here the sign of the gauge field 



strength A is not important, since the physical quantities de- 
pend only on the parameter A? as we will show in this paper 
In this paper, we set /i > without loss of generality. 

We now turn on a short-range attractive interaction in the 
spin-singlet channel. For homogeneous Fermi gases, we de- 
fine the Fermi momentum k^ through the fermion density 
n - N/V = k^l^iir), and the Fermi energy is ep = k^jl. 
In the dilute limit ^f'"o ^ 1 ('"o- effective range of the in- 
teraction), the interaction Hamiltonian can be modeled by a 
contact interaction. The total Hamiltonian of the system can 
be written as 



fl'ViA"'"(r)C7/o+Ho)^^(r) 



U \ dh ilj'\{r)il,\{r)ilf^{r)ilj^{r), (3) 



where 'Ko = p~/2 - yU is the free single-particle Hamiltonian 
with ji being the chemical potential, 'Kso = /Ir ■ p is the SOC 
term, and U > denotes the attractive s-wave interaction be- 
tween unlike spins. For the validity of such a contact interac- 
tion, another dilute condition /Iro <K 1 should be satisfied ll29tl . 
We also note that the Galilean invariance of the Hamiltonian 
in the absence of SOC {A - 0) is broken by the SOC term. 
However, as it will be shown, the Galilean invariance can be 
recovered at the boson (rashbon) level for large A. 

In the functional path integral formaUsm, the partition func- 
tion of the system is 



where 



Z = J !Di//Di{fe\p{-S[tfr, tfr]}, 

rP r rfi 
Jo J Jo 



(4) 



•A). (5) 



Here = 1/T and H{t//, iff) is obtained by replacing the field 
operators and iff with the Grassmann variables and if/, re- 
spectively. To decouple the interaction term we introduce the 
auxiliary complex pairing field <i>{x) = -U4'i{x)4'-\{x) {x - 
(r, r)] and apply the Hubbard-Stratonovich transformation. 
Using the 4-component Nambu-Gor'kov spinor '^{x) - 
{'^■\, tf/i, iAt, lAi]^' we express the partition function as 

Z ^ J OTOTSOBO* exp I - J" dx\Q>(x)\^ 

+ ^ J dx J dx''i'{x)G-\x,x')'¥ix')^, (6) 

where the inverse single-particle Green function G '(x, x') is 
given by 



G ' = 



I _, -dr-'Ho- 'H,o iTy^(x) 



, \6(x-x'). (7) 



Integrating out the fermion fields, we obtain Z - 
J £)<1)£)<1)* exp { - .Seff[0,<l)*]}, where flie effective action 
reads 

SM dxmxt - ^Trln[G-'(x, x')]. (8) 
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III. TWO-BODY PROBLEM 

In this section, we study the two-body problem at vanish- 
ing density. We will determine the binding energy and effec- 
tive mass of difermion bound state formed in the non-Abelian 
gauge field. The systematic way to study the two-body prob- 
lem in presence of a nonzero spin-orbit coupling A is to con- 
sider the Green function r{Q) of the fermion pairs, where 
Q - {iv„,q) with v„ - 2nnT (n integer) being the bosonic 
Matsubara frequency. For zero density, we need to consider 
the case = 0. In the functional path integral formalism, 
r '(G) can be obtained from its coordinate representation de- 
fined as 



1 52^eff[<I>,<I>*] 



(9) 



ygy 5<t)*(jc)(50(x') 

For O = 0, the single-particle Green function G{K) reduces to 
its non-interacting form 



g^K) 
g-(K)]' 



(10) 



where K - (iw„, k) with cl>„ = (2« + l)nT being the fermionic 
Matsubara frequency. Here the elements g±{K) read 



gAK) = 
g-(K) = 



1 



iaj„ - f k - ^so 
1 



(11) 



ia)„ +^k- ^so 

Here = ek - A' with ek = ^so = '^(T.v^.r + Tyky + r-kj 
and ^*(, = AiTfkx - Tyky + r-kj. The inverse in g±{K) can be 
worked out and we obtain 



gAK) = 
g-(K) = 



ia)„ - ^k + ^so 

(/w„-^k)2-i2k2' 

(/w„+^k)2-^2k2- 



(12) 



The single-particle excitation spectrum therefore has two 
branches, = ^k ± A\k\, due to the spin-orbit coupling. 

Using the free fermion propagators g+{K), F '(0 can be 
expressed as 



(13) 



Completing the Matsubara frequency sum and making the 
analytical continuation /v„ oj + iO^, the real part of 
F"'(w -H iO^, q) takes the form 

FR'(w,q) = ReF-'(w-H/0+,q) 

1 1 V ^ ~ /(^k+q/2^ ~ /(^k-q/2^ /, ^ \ 
(l+«r^kq), 



f/ 4 ^ , 

o-,y=± k 



^k+q/2 ^ ^k-q/2 



(14) 



where f(E) = 1 / {e^^ -)- 1 ) is the Fermi-Dirac distribution func- 
tion, and Tkq is defined as 



kq 



(k + q/2)-(k-q/2) 
|k + q/2||k-q/2| 



(15) 



We use the notations YiK - T Ym Sk ^"d 2k = / d^^Kln)^ 
throughout this paper. Note that F~'(Q) takes the form similar 
to that of the relativistic systems lUHl . due to the fact that 'H^o 
behaves like a Dirac Hamiltonian. 

The integral over the fermion momentum k is divergent and 
the contact coupling U needs to be regularized. For a short 
range interaction potential with its s-wave scattering length a 
it is natural to regularize U by means of the two-body problem 
in the absence of SOC. We have 



1 

U 



1 

Ana^ 



y-- 

^ 2ek 



(16) 



For the pure two-body problem at vanishing density and 
temperature, we discard the Fermi-Dirac distribution function. 
The energy-momentum dispersion coq of the pair excitation is 
defined as the solution oj + l/j. = aiq of the two-body equa- 
tion Fj^' (a>, q) = . After some manipulations, the two-body 
equation becomes 



z 



1 

k2 



^kq 



£} -4^2k2 
kq 



4/1^ k^q- sin- 



1 



Anas 



(17) 



Here (p is the angle between k and q, and fikq = fk+q/2 + 

ek-q/2 - Wq = k^ H- q-/4 - Wq. 



A. Bound State and Binding Energy 

We are interested in whether there exist difermion bound 
state in the presence of SOC. For this purpose, we first con- 
sider zero center-of-mass momentum q and determine the en- 
ergy regime where the imaginary part of r^^(a) + iO^, q = 0) 
vanishes. We have 



ImF-'(w-H/0+,q = 0) 



— -y r 



k^dkdik^ + laAk ■ 



■2fi). (18) 



Therefore, a bound state exists if the equation Fj^' (o), q = 0) = 
has a solution in the regime -oo < u + Iji < -A?. 

The binding energy in the presence of nonzero SOC 
is determined by the solution of w -i- 2/z - -Eb for the 
equation Fj^'(w, q = 0) = 0. From the imaginary part of 
F"'(dLi + ;0^,q - 0), the binding energy E^ must be larger 
than a threshold isth = 't^- The equation determining reads 



r 

Jo 



k^dk 



k^ + Ef, 



1 



k^- (k^ + £b)2 - 4^2^2 



- f . (19) 
2a 



Completing the integrals analytically, we obtain a simple al- 
gebraic equation for E^, 



Eb - 2A- 1 



(20) 



We find that, for arbitrary scattering length a,, there always 
exists a solution Eb > A^ . Therefore, the difermion bound 
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state can form in the presence of SOC even for a, < where 
no bound state exists in the absence of SOC. 

The solution of Eq. ( l20l i can be analytically expressed as 



(21) 



Therefore, the quantity Eb / A~ depends only on the dimension- 
less parameter k = l/iAus). We have 



where the function J'(k) is defined as 



(22) 



(23) 



We are interested in the case Aos — > oo or /c = 0. This happens 
when a , — > oo (unitary point of the Feshbach resonance) for 
fixed /I or /I — > oo for fixed Qs- In this case, we have J' - 2 
and a very simple result 



EB(Aas 



oo) = 2r 



(24) 



In general, the numerical result for the quantity E^l A? - 1 
J^{k) - 1 is shown in Fig. [1] 




FIG. 1: The quantity E^/A^ - 1 as a function of the dimensionless 
parameter ^= l/(/ia,). 



Defining a new variable x - k/A, this equation becomes 

2 (x-+jf+4-x^ 



_ 2m\ r 
'"b/ Jo 





dxx 



[(x2 + J")2 - 4x2]2 



(x2+J")[(x2+ J)2-4x2]2- 

Completing the integrals analytically, we obtain 



2m 

OTb 



4/' J'-l 
3l 



3/2 



2 



(26) 



(27) 



The eff'ective mass therefore depends only on the combined 
parameter A- = l/{Aas). We have 



2m 
OTb 



7 



3 + 4 + K ^K^ +4 

4 / 2 

3 \ 2 _,_ 4 ^ V/f2 + 4 



3/2 



(28) 



The numerical result for ;7Zb/2ot is shown in Fig. |2] We find 
analytically that m^ — > 2ot in the limit a: — > +oo and otb — » 6m 
in the limit a: — > -oo. For the case /ia, — > oo or a = 0, the 
effective mass reads 



OTB(/lfl.5 -» oo) 3(4 + V2) 



2m 



14 



= 1.16. 



(29) 




B. Molecule EfTective Mass 

For small nonzero center-of-mass momentum q, the solu- 
tion for dLiq can be written as = -E^ + q'/{2mB), where 
Otb is referred to as the effective mass of the bound state. Sub- 
stituting this dispersion into the equation rj^'(w, q) = and 
expanding the equation to the order 0(q'), we obtain 



k^dk- 



[(/t2 + £b)2 - 4A^k^Y 
U^k^ 

(fc2+£B)[(fc2+£g)2_4^2^2]2- 



(25) 



FIG. 2: The molecule effective mass ;;ib (divided by 2m) as a func- 
tions of the dimensionless parameter k = 1/ (Aas). 



In summary, the difermion bound state forms in the pres- 
ence of SOC for arbitrary value of s-wave scattering length 
a.5. The bound state possesses a non-trivial binding energy 
^B and a non-trivial effective mass otb > 2ot. Such type of 
bound state are referred to as rashbon in the previous litera- 
tures. Due to the formation of bound state at the BCS side of 
the resonance (a^ < 0) and the enhancement of binding en- 
ergy in the presence of SOC, we expect that there will be a 
crossover from the BCS superfluid state to the Bose-Einstein 
condensation of rashbons if the spin-orbit coupling A can be 
tuned from small to large values. 
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IV. SUPERFLUID GROUND STATE: MEAN FIELD 
THEORY 

For the many-body problem, we first consider the properties 
of the superfluid ground state {T = 0) in the self-consistent 
mean-field theory. In the superfluid ground state, the pairing 
field <i>{x) acquires a nonzero expectation value (<l>(x)) = A, 
which serves as the order parameter of the superfluidity. With- 
out loss of generality, we set A to be real. Then we can express 
the pairing field as = A H- where <p(x) is the fluctu- 
ation around the mean field. The effective action »Seif [O, O*] 
can be expanded in powers of the fluctuation. 



5efF[1>,1>*] 



(30) 



where »S|,^(A) = >Seif[A, A] is the saddle-point or mean-field 
effective action with the superfluid order parameter deter- 
mined by the saddle point condition dSf^/dA = 0. 

In the mean-field approximation, the grand potential Q. = 
»SefF[A, A]/(fiV) can be expressed as 



1 1 



(31) 



where the inverse fermion Green function reads 



Since the integral is convergent, we can use the trick k^+2Ak = 
{k+A)^-A^ and convert the integration variables to k+A. Then 
we obtain 

C" dk T T I A? . . \ 



where 



= £k - fi, Ek 



(38) 



with fi - fi + A? 12. 

Using the above expression for Eq, the gap A and the 
chemical potential can be determined by dEc/dA = and 
dEcldji - -n, i.e.. 



f 

Jo 

f 

Jo 



dkik^ + A^) 



dk(k^ + A^) 



^(Q-/i)-+A2 



= 0, 



(39) 



(32) 



Using the formula for block matrix, we first work out the de- 
terminate and obtain 

det^-i(/^^«, k) = [{iaj„f - (£+)2] [{ico„)- - (E^)'] , (33) 



where E^ - ^/(fir^Tlikjj^H^'A^ are quasiparticle excitation 
spectra. Then completing the Matsubara frequency sum we 
obtain 

" = ^ + Xi[^''-nv(£+)-nv(£k)], (34) 



where "WiE) = E/2 + rin(l + e''^''^). Note that the term 
Sk = 5 I]k(^k ^k ^ added to recover the correct ground 
state energy for the normal state (A = 0). 



We notice that the above expressions for the gap and num- 
ber equations can be analytically evaluated using the elliptic 
functions, like the analytical treatment for the gap and number 
equations in the absence of SOC ll30ll . 



B. Fermion Green Function 

The explicit form of the fermion Green function Q{i<x)„, k) 
can be evaluated using the formula for block matrix. In the 
Nambu-Gor'kov space, it takes the form 



^(iw„,k) 



^ii(/w„,k) ^i2(/w„,k) 



(40) 



The matrix elements can be expressed as 



A. Ground State Energy 

At zero temperature, the ground state energy Eq = Q,(T - 
0) is Ec = A^/U + (1/2) 2k(2fk - £k - -^k)- Using the fact 
that the binding energy satisfies the equation 



U 2 ^ 



■DO a 



k^dk 



1 



27r2 k^ + 2a Ak + Eb ' 



(35) 



guim,, k) = y[ii(/w„, k)/ + Sii(/w,„ k)M, 
Qiiiioj,,, k) = y[22('w„, k)/ + S22(/w„, k)M*, 
gniiojn, k) = -iTy [^i2(/w„, k)/ -t- Snim,, k)M*] , 

giiimu k) = iTy [^21 (icj„, k)/ + S21 (m„ k)M] , (41) 

where / is the identity operator in the spin space and the oper- 
ators M and M* are defined as 



we can express the ground-state energy in terms of as 



^g-2ZjJo '^\k^+2aAk + Es ' ^ 



(36) 



M 



M 



^ Txkx — Tyky + T~k- 



(42) 
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The explicit form of the quantities Jljj and Sij are given by 



C. Gap and Chemical Potential 



^ii(/w,„k) = i 



J?l22(!W«,k) = - 
^12(!W«,k) = i 



ioj,, + 



^2i('w„,k) = J?[i2(/w„,k) 



(43) 



and 



Sii(iw„,k) = - 



(!W„)2 - (/W„)2 - (£^)2 



S22(iW„,k) 



Si2(!w„,k) 



1 



S2l(!W„,k) = -Si2(/W„,k). 



(44) 



Using the matrix elements of the Green function, we can 
calculate various quantities. First, the momentum distribu- 
tions «|(k) and «4,(k) for the two spin components can be eval- 
uated as 



= -y 



^ii(iw,„k) + -^Sii(idLt,„k) 



= -y 



J?[ii(/w,„k) - -^Sii(/w,„k) 



(45) 



Second, the singlet and triplet pairing amplitudes can be ex- 
pressed as 

'^n(k) = (lAkTiA-ki) 



= -y 



-yi2i(iw„,k) + -^S2i(/w,„k) 
k 



0iT(k) = <lAkilA-kT> 



- -y 
p^ 



yi2i(/w,„k) + -^S2i(/w„,k) 
k 



(^ll(k) = {lAkTlA-kT) 

^.t - ifcy 1 



^ J]S2i(/.;„,k), 
<^ii(k) s (lAkiiA-ki) 

Third, the gap equation for A can be expressed as 
A = -C/i^^y[i2(/o.„,k). 



(46) 



(47) 



Using the ground state energy Eq, the original forms of the 
gap and number equations at T = are 



U 22(2^+ +2£,:)' 



(48) 



The pairing gap A and the chemical potential ^ can be numer- 
ically solved for given values of 1 /{kpOs) and A/kf. From now 
on, we denote the saddle point solution for the gap at zero 
temperature as Aq. We also notice the relation 



1 

Aa,, 



kpGs \kf 



(49) 



(A) Analytical Results for Large SOC. We first obtain the 
analytical solution at large SOC, A/kp » 1. For large SOC, 
we expect < and Aq <sc |yu|. Therefore, we can expand the 
equations in powers of Ao/|/i| and keep only the leading order 
terms. The gap equation becomes 



k^dk 



1 



27r2 + 2a Ak - l/j 



\2 N 



Comparing with the two-body problem, we obtain 



' 2 ' 



(50) 



(51) 



Substituting this into the number equation, we obtain 

A2 



_ ^ r 

^2 Jo 



k'dk 



1 



1 



+ 



k ' ^^A- • 



k dk 

[(F -t- £b)2 - 4^2^]^ 



(52) 



We notice that this integral also appears in Eq. (25). Complet- 
ing the integral analytically, we obtain 



^ AnAn 



(J- -1)3/2 



J 

AA [2eF(J--l)F 



Therefore we have 
Ao 



16 (J - 1)3/2 ^ 



J h 

In the limit Aa^ — > 00, we have 3^-2 and therefore 

2^(2eF)3/2 



ts.Q{Aas — > 00) ^ 2TiAn 



In 



(53) 



(54) 



(55) 
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FIG. 3: The pairing gap Ao (divided by ep) as a function of A/kp. The 
red dashed hne shows the analytical result (54) or (60). 
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FIG. 4: The chemical potential yu (divided by ep) as a function of 
A/kf. The red dashed line shows the analytical result ^ ^ for 
large SOC. 
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It can be written as another interesting form 



8 A 

37T kp 



(56) 



Therefore, for very large SOC, the gap Aq increases as Aq ~ 

VI. 

Beyond the leading order, we can write the chemical poten- 
tial jj as 



(57) 



where fi^ - Ifi+E^ <g: Eg is referred to as the effective chem- 
ical potential for bosons (rashbons). We will give an explicit 
expression for /iB in Section V. 

(B) Numerical Results. The gap and number equations 
(39) and (48) are equivalent. For numerical calculations, it is 
convenient to employ Eq. (39). If we define the following 
dimensionless quantities 



1 



Xl = — , 



X2 



Ao 



— , (58) 



the gap and number equations can be written as the following 
dimensionless form 







Jo 



1 



1 



4 



(z'-g 



dziz^'+gl) 



■Xi)2 
_ ^2 



-Xl 



^(z2-^2-Xi)2 + x2 



= ^.(59) 



The integrals in the above equations can be analytically eval- 
uated using elliptic functions l3Qil . For given values of gi and 
g2, these two equations determine xi and X2. 

The numerical results are shown in Fig. [3]and Fig. |4] The 
red dashed lines correspond to the analytical results for large 
SOC, 



Xl = -giJigi/gi), 



Xl = 



' i6gi mg2igi)-i?i^ 

37T J{g2lgl) 



(60) 



We find that the pairing gap generally increases with increased 
/l/fep. The numerical results become in good agreement with 
the analytical results when /l/fep > 1. Therefore, the system 
enters the rashbon BEC regime at Ajk-p ~ 1 . For large positive 
value of l/ikpOs), the analytical results are in good agreement 
with the numerical results even for small values of A/kp. For 
very large A, we find the numerical results fit very well with 
the following scaling behavior 



Ao 



I A 
37r V fcp ' 



^ -2, 

ep \kF 



A 



(61) 



for both negative and positive values of l/{kpas). 



9=0 



9=Ji/3 



0.5 1 1.5 2 2.5 3 3.5 4 




FIG. 5: The fermion momentum distributions nf{k) and for 
various values of the polar angle 6. We set l/(fcpa,) = -1 and A /hp = 
1 in this calculation. 



D. Fermion Momentum Distribution 



From the matrix elements of the fermion Green function 
0{ia)n,k), we can obtain the momentum distributions «-r(k) 
and «x(k) for the two spin components. The density of each 
component reads «o- = Y^k n(T(k).We find that even though the 
density of the two components are the same, «-f = n^, their 
distributions in the momentum space are different. At zero 



9 



temperature, their explicit expressions are given by 

«KM, = iZ(.-|)^T?''(-|). 

"*« = iz(.-|)-T?"(-|).<-> 

where 6 is the polar angle in the momentum space. We find 
that n|(k) - «x(k) only for G - n/2. We have n-f(k) < «j,(k) 
for < < 7T/2 and «|(k) > nx(k) for 7t/2<6< n. 

In general, with increased SOC, the distribution broadens, 
which indicates a BCS-BEC crossover. A numerical exam- 
ple for \l{k-pas) - —\ and Ajk-p = 1 is shown in Fig. |5] 
The new feature here is that the distributions generally display 
non-monotonous behavior due to the SOC effect. We note that 
the peaks in the distributions are just located ?Xk - A. 




0.2 



0.5 1 1.5 2 2.5 3 

FIG. 6: The condensate fraction 2No/N as a function of A/kp for 
various values of l/(A:Fa.,)- 



E. Condensate Density 



According to Leggett's definition H3Ifl . the condensate num- 
ber of fermion pairs is given by 



Yj // «?'ri£/^V2K<A<r,(ri)«A^,(r2))|-. (63) 



0-l.CT2=T,i 



For systems with only singlet pairing, this recovers the usual 
result A^o = // d^rid^r2\{tfr^(ri)if,i(r2)}\^ (H. Converting 
this to the momentum space, we find that the condensate den- 
sity «() = No/V is a sum of all absolute squares of the pairing 
amplitudes. 



no 



i ^ [|<^n(k)P + \<Pim- + l0TT(k)P + l<^u(k)Pl 



Tj \ ^2'^2iOw«,k) + i J]S2i(/w„,k) 



y6- 



(64) 



Completing the Matsubara frequency summation and tak- 
ing the zero temperature limit, we obtain the explicit expres- 
sion for T - Q, 



no 



Jo 



167T- 



k^dk 



1 



1 



e + A^ 



(Q-Ai)' + A, 



(65) 







Generally, we can show that no < n/2. For large SOC and/or 
attraction, we have Aq <*c Using the number equation (39) 
or (48) and expanding all terms in powers of Ao/|ju|, we find 
that 



no---0 



(66) 



Therefore, the condensate fraction 2No/N approaches unity at 
large SOC and/or attraction, as we expected. 



In general, the condensate fraction 2No/N can be expressed 



as 



2No 
N 



3x1 



f 

Jo 



dz- 



(67) 



It can be numerically obtained using the solutions of x\ and 
X2 from the gap and number equations. The numerical results 
are shown in Fig. |6] We find that, even for negative values 
of 1 / (kpas), the condensate fraction approaches unity around 
A/kp ~ 2. This is consistent with the observation from the so- 
lutions of the gap and number equations that the system enters 
the rashbon BBC regime at A/kp ^ 1 for negative and small 
positive values of I /{kpas). 



F. Superfluid Density 

To evaluate the superfluid density «j, we can employ the 
standard definition ll33il34ll . When the superfluid moves with 
a uniform velocity = {v^, u,,, i;.), the superfluid order pa- 
rameter transforms like ^ Oe^"!' and O* cD'e ^'q, 
where = mVs (m = 1 in our units). The superfluid density 
« 5 is defined as the response of the thermodynamic potential 
Q. to an infinitesimal velocity velocity v,, i.e.. 



Q(q,) = Q(0) + i«,q2 + O(q^). 



(68) 



The thermodynamic potential in the presence of a velocity Vs 
can be evaluated by a gauge transformation for the fermion 
field i// i/rg-'i' ■•. We have 



f^(q..) lndet^;'(/o.„, k). 



(69) 



II k 



where the inverse fermion Green function in the presence of 
Vs reads 



g-'(icu„,k) = g-\iio„,k) - Z(q,). 



(70) 
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Here the velocity-dependent part S(qi) includes three parts, 
I.{qs) = 2:i(q,) + S2(q.O + S3(q.,), where 

1 7 

Si(q.) = 2^sO-^' 
22(q.) = k ■ q,l, 

5^3(q.5) = /1(t^4.vO"3 + Tyqyl + T.q-a-3). (71) 

Here cr,- (/ = 1, 2, 3) and 1 are the Pauli matrices and the iden- 
tity matrix in the Nambu-Gor'kov space, respectively. We 
note that the term S3(q.O is purely due to the presence of SOC. 

(A) Derivation of the Superfluid Density. The superfluid 
density can be obtained by the method of derivative expan- 
sion for Q(q.5), i.e., 

Q(q,) = a(0) + ^ V -Tr [^S(q,)]" . (72) 
2 ^ n 



We find that there are four types of nonzero contributions at 
the order (9(qj): 



Qi ~ Tr(^Ii), Q2 ~ Tr(^22^22), 

n^~Tr{gz^g^i), ii,~Tx{g^2S^i)- (73) 



Since the superfluid state is isotropic, the superfluid density 
tensor should also be isotropic. We have carefully checked 
that all anisotropic terms vanish exactly. Completing the trace 
in the Nambu-Gor'kov and spin spaces, we finally obtain the 
following expressions for the four types of contributions; 



n k 



n k 



2 P 



11^11 



- ^22^22 + 2^2lS2l) . 



(74) 



H k 



Note that the first contribution is just from the total particle density n, Q.i - ^nq^. Collecting all terms, the superfluid density «j 
is given by 



= n+ ^XiZ y + + ^22 + ^li + 2^21 + 2S^i) + ^ (^uSn - ^22^22 + 2^21^21) 



n k 



(^tl + -5^22 + 2^21) - -T (^11 + ^22 + 2S^i) 



(75) 



Completing the Matsubara frequency sum, we obtain the finite-temperature expression 



r°° k^dk 

_ /( f" kdk 
" 3 Jo 2^^ 



(/fc-i-A)2 1 



(A:-/l)2 1 



6 2r^^^,2 



6 2r^^^,2 



A2\1-2/(£,+) 



A2\1-2/(£-) 



(76) 



We have checked that this expression is consistent with the re- 
sult for ordinary fermionic superfluids in the absence of SOC 
miH]. Also, setting A = 0, we find that ns(A = 0) vanishes 
exactly. 



We are interested in the zero temperature case. At zero tem- 
perature, the superfluid density reduces to 



(77) 
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where n,i is given by 

A r 



kdk 



(78) 



We notice that n i vanishes in the absence of SOC and we re- 
cover the usual result « , = « at T = for ordinary fermionic 
superfluids lf33l[34ll . However, for nonzero SOC, is always 
positive and we have risiA + Q) < n. Therefore, the SOC leads 
to suppression of the superfluid density. 

(B) Analytical Result for Large SOC. To understand this 
interesting phenomenon, we first take a look at the large SOC 
limit. In this case we have yu ^ -Eb/2 and A <sc Therefore, 
we can expand the expression in powers of A/|yu| and keep only 
the leading order terms. Doing so, we obtain 



A2 

_I2. 

87r2 

A," 



k^dk 



1 



1 



Jo 



k' ^''k' 
■2 I i7_\2 , /I i27,2 



k dk 



(79) 



and 



67r2 Jo U, 



1 



A^4 
7r2 3 



k^dk- 



1 

" 2 



(^2 + Eb) [{k^ + £b)2 - 4A^k'-] 



(80) 



Comparing the above results with the equation for the 
molecule effective mass, we find that n^/n ^ 1 - 2m /ms- 
Therefore, at large SOC, the superfluid density is suppressed 
by a factor 2m/ms, i.e.. 



2m 
ris - — n. 
ms 

For A ^ oo, using the result for 2m /m^ at k 
the ratio n jn approaches a universal value. 



(81) 



0, we find that 



-Wkp 

n 



oo) 



14 



3(4 + V2) 



= 0.862. 



(82) 



To further understand this result, we consider the eff'ec- 
tive action for the phase field G{x). To this end, we write 
the order parameter as <l)(x) - A(x)e'**'^*. In the static limit, 
we can obtain the effective Hamiltonian for the phase field, 
i/gff - (JJ2) J d^r[Vd{r)]^, where the superfluid phase stiff- 
ness Js is related to the superfluid density by 7, = n j{4-m). 
Therefore, at large SOC, we have 



J 2m n 
' ms 4m mB ' 



(83) 



where hb = n/2 is the density of bosons (rashbons). This 
means that, at large SOC, the superfluid phase stiffness self- 
consistently recovers that for a rashbon gas with a non-trivial 
effective mass m^. We emphasize that this interesting result 
was first observed by us in 2D Fermi gases with Rashba spin- 
orbit coupling ll25ll . 



1 . 
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FIG. 7: The superfluid density n, (divided by n) as a function of A/kp 
for various values of l/(fcFaj). The red dashed line shows the ana- 
lytical result 2m I m^. The dot-dashed line corresponds to the value 
0.862 for /Iflj — » oo. 
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This result also indicates that the Galilean invariance, 
which is explicitly broken in the original fermion Hamilto- 
nian, can be effectively restored at the boson (rashbon) level 
at large SOC. This is due to the fact that at large SOC the sys- 
tem becomes a weakly interacting Bose-Einstein condensate 
of non-relativistic rashbons which have a non-trivial effective 
mass niB- We will show this conclusion explicitly in the next 
section by deriving the Gross-Pitaevskii free energy for the 
dilute rashbon condensate at large SOC. 

(C) Numerical Results. The superfluid density at zero tem- 
perature can be expressed in terms of the dimensionless pa- 
rameters as 



z + 2agiz - xi + ^ 



(84) 



2agiz - xi)2 -H xl 



It can be numerically obtained using the solutions of xi and 
X2 from the gap and number equations. 

The numerical results for n jn as a function of A/kp for 
different values of l/(kpas) are shown in Fig. |7] For nega- 
tive values or small positive values of l/{kpas), the numerical 
result becomes in good agreement with the analytical result 
^ Im/niB when A/kp > 1, which is consistent with the 
observation that the system enters the rashbon BBC regime at 
A/kp ~ 1. For large positive values of l/ikpOs) (in fact even 
for l/(kpas) = 1), the numerical results are always in good 
agreement with the analytical result for all values of A/kp. 

For both negative and positive values of 1 / (kpOs), we find 
that n/n, approaches a universal value 0.862 when A/kp — » oo, 
as indicated from the analytical observation. 



A. Calculation of the Ginzburg-Landau Coefficients 

The coefficients c and d of the potential terms can be ob- 
tained from the mean field thermodynamic potential Qq = 
(r/y).Seif[A' , A] which can be evaluated as 



1^1 _ y 
Ana. ^ 



2ek 



(86) 



We have 



c - 



5-Q 



A=0 



5(|A|2)2 



(87) 



A=0 



After a simple algebra, the coefficients a and /3 can be evalu- 
ated as 



c = 



1 

4^ 



' -2jj -2A^ r 

1 -2// -I- 2A^ 
16^ (-2fi - A^y 



/2- 



(88) 



From the expression of c, we find that a quantum phase 
transition from vacuum to Bose condensation takes place at 
- -Eb/2. Thus near the phase transition, c can be simpli- 
fied as 



1 



8^(£b-^2) 



3/2 



(89) 



V. BOSE-EINSTEIN CONDENSATION OF WEAKLY 
INTERACTING RASHBONS 

As we have shown in the last section, the superfluid state in 
the large SOC limit is a Bose-Einstein condensation of rash- 
bons. We are interested in the interactions among the rash- 
bons. In this section, we will derive the Gross-Pitaevskii free 
energy for a dilute rashbon condensate, which allow us to 
extract the rashbon-rashbon scattering length. Another goal 
of this section is to show that the Galilean invariance, which 
is explicitly broken in the original fermion Hamiltonian, can 
be effectively recovered at the boson (rashbon) level at large 
SOC. 

To this end, we consider the mean field theory where the 
auxiliary boson field <l)(x) is replaced by its expectation value 
{(i>{x)) - A{x). In the large SOC limit A ^ oo, the fermion 
chemical potential fi approaches -Ep,/2. Since the pairing gap 
|A| <sc |yu|, we can expand the effective action in powers of 
|A| (as well as in powers of its space-time derivatives), which 
results in a Ginzburg-Landau free energy functional 



VodMx)] 



= ^ dx AHx)[a^-bV^'^A{x) 



+ c|A(x)|2 + -d\Mx)t 



(85) 



where jUb = 2/i -I- Eb ^ £b is the boson chemical potential. 
Further, setting - —Ep,/2, d can be reduced to 



1 Eb + 2A- 



(90) 



The coefficients a and b of the kinetic terms can be obtained 
from the inverse boson propagator with A = 0. It can 

be evaluated as 



1 1 



1 + ajT ,q 



^ ^± ^ ^k+q/2 + ^k-a/2 



(91) 



In the large SOC limit, the coefficients a and h can be obtained 
by the small momentum expansion for Tr^{Q). We have 



2)"'(e)--fl(iV„+//B-^|, 



(92) 



where mp, is the rashbon effective mass determined by (28), 
and a is given by 



1 



87r (£b -^2)3/2 ■ 
We observe the relation c - £)"'(0) = -a^B- 



(93) 



13 



B. Gross-Pitaevskii Free Energy 

According to the above results for the Ginzburg-Landau co- 
efficients, if we define the new condensate wave function \\iix) 
by 



(94) 



the Ginzgurg-Landau free energy can be reduced to the Gross- 
Pitaevskii free energy of a dilute Bose gas. 



Vgp[v(^ 



x)] ^ Jdx V'(Jc)^^ - ^j^' 

2 1 47rflBB , , .,4 

/^b|v|/(-^)I + i¥(-^-)r 

2 ms 



(95) 



where abb is the boson-boson scattering length. Its explicit 
expression is 



Eb + 2A^ r- 

Abb = wJb -5 V-Eb - 



(96) 



Note that m = 1 in our units. For A - Q and > 0, us- 
ing the result otb - 2 and ^b = l/a,, we recover the well- 
known result Abb = 2a s iH). One remark here is that this re- 
sult is the mean field result which is not exact. In the absence 
of SOC, exact four-body calculation shows that obb - 0.60 5 
ll35ll . Therefore, it is interesting to explore the exact rashbon- 
rashbon scattering length in the future studies. Another theo- 
retical framework to obtain more exact abb is to include the 
Gaussian fluctuations ll36ll . 

The Gross-Pitaevskii free energy explicitly shows that the 
Galilean invariance, which is explicitly broken in the origi- 
nal fermion Hamiltonian, can be effectively recovered at the 
boson (rashbon) level at large SOC. 



C. Rashbon-Rashbon Scattering Length 

Using the expressions for the binding energy Es and the 
effective mass niB, we obtain 



Abb = - 

A 



1 2(J + 2)^/J^ 



r- 



7 _ 4 (Jzl\^'^ _ 2. 
3 i \ J I 3 



(97) 



We find that the quantity /Iobb depends only on the dimen- 
sionless parameter k - \j{Aa^. For the case Aus ^ 00 or 
/<■ = 0, we have 3^-2 and therefore 



flBB(/lfl.! 00) 



1 3(4 + V2) 2.32 



A 



(98) 



The numerical result for the scattering length abb is shown in 
Fig. [8] We find that the quantity Aa^^ has a maximum near 
the point /c = 0, at = -2. 11 . 




FIG. 8: The molecule scattering length obb (divided by IjX) as a 
function of the dimensionless parameter k = 



D. Rashbon Chemical Potential 

For a uniform system, the expectation value of the con- 
densate \\i{x) should be determined by minimizing the Gross- 
Pitaevskii free energy. We find that the minimum is given by 



^0 



(99) 



where go - Ana-R-RlniB- The total density of the bosons is 
„B = «/2 = IvoP 
potential can be given by 



oAq. Therefore, the boson chemical 



Zttobb 
jUB = n. 



(100) 



For the case Aa^ 
obtain 



00, using the result for and a^B, we 



fiBiAa,, -> 00) 



Inn 



(101) 



VI. GAUSSIAN FLUCTUATION AND COLLECTIVE 
EXCITATIONS 

To study the collective excitations, we consider the fluctu- 
ations around the mean field. Making the field shift <I)(jc) — > 
Ao + (f>{x), we can expand the effective action S^s in powers 
of the fluctuations. The zeroth order term S'"^^ is just the mean 
field result, and the linear terms vanish automatically guaran- 
teed by the saddle point condition for Aq. The quadratic terms, 
corresponding to Gaussian fluctuations, can be evaluated as 

'5lff*[<^,0"'"] = jXl-^'^e) <^(-G))M(0(^fJ® Al02) 



where the inverse boson propagator M takes the form 
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with the relations Ml i(0 - M22(-0 andMi2(0 = M2i(0. At zero temperature, the explicit form of M(0 can be evalu- 
The matrix elements of M.(Q) can be expressed in terms of the ated as 
fermion propagator ^(A"). We have 

Mii(0 = - ^ y Tr[gn(K + 0Tv^22WTy] , 
Mi2(0 = ^ 2 + Q)rrgn{K)T,] . (104) 

I 



(«Lq/2) («k-q/2)" (^'k-q/2) 



k+q/2 k-q/2 



^ ^k+q/2 ^ k-q/2 



(l+ayTkq) 



(105) 



and 



Mi2(0 



a,y-± k 



;/"■ v" 1/ y''' 
"k+q/2 ''^k+q/2"k-q/2 ' k-q/2 



y"" y'' 
"k+q/2 ''^k+q/2"k-q/2 '^k-q/2 



+ ^k+q/2 



k-q/2 



^k+q/2 



k-q/2 



(1 + ffyTkq) . 



(106) 



Here the BCS distribution functions are defined as (v^)^ - 
(1 - ^'l/E'^)/! and (m^O^ = (1 + ^'^/E'^)/2. In the absence of 
SOC, A = 0, the expressions for Mii(2) and Mi2(0 recover 
the results obtained in jst]. 



A. Bogoliubov Excitation in the Raslibon Condensate 

At large SOC and/or attraction, the superfluid state is a 
Bose-Einstein condensation of weakly interacting Rose gas. 
Thus we expect that the low energy collective excitation in 
this case recover the well known Bogoliubov excitation spec- 
trum in a weakly interacting Rose condensate iIjtIi . In this 
part, we will give an explicit proof for this. 

In the large SOC and/or strong coupling limit, the chemical 
potential reads /i ^ -Eb/2 and we have Ao «: In this case, 
we can expand the matrix elements of M in powers of Ao/|yu| 
and keep only the leading-order terms. Following this spirit, 
we obtain 



Mn(Q)^r>-\Q)+XAl, 
Mi2(0 - YAl 

where the coefficients X and Y are given by 



1 



1 



2d. 



(107) 



(108) 



boson propagator M{Q) can be well approximated by 



Mii(0 ^ a \-iv„ + ^ - + 2golVol^J 

M22(0 - a ^iv„ + ^ - yUB + 2gol¥ol"j 
Mi2(0 = M2i(0-fl§olVol', 



(110) 



where go = Ajia^^jm^ and |\|/o|" - fJ^a/go is the minimum of 
the Gross-Pitaevskii free energy (corresponding to the saddle 
point Ao of the effective potential). From the Gross-Pitaevskii 
free energy, the boson density reads «b = n/2 - |\|/oP- Utiliz- 
ing these results, we obtain 



M(0 - a 



+ 2^ + ^o"B 



gons 



, (111) 



Taking the analytical continuation iv„ — » oj+iO^, the disper- 
sion a> = a;(q) of the collective mode is obtained by solving 
the equation 



detM[q,w(q)] = 0. 



(112) 



Therefore, the Goldstone mode takes a dispersion relation 
given by 



Further, taking the small momentum expansion for D '(2), 
we obtain 



Mn(0 - -a\iv„ +i^B- ^j + 2JAg. 



(109) 



Therefore, in the large SOC and/or strong coupling limit, the 



w(q) 



2otb \2mB niB /' 



(113) 



This is just the Rogoliubov excitation spectrum in a dilute 
Rose condensate where the bosons possess a mass and a 
two-body scattering length abb- 
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B. Collective Modes in the BCS-BEC Crossover 

The dispersions of the collective modes are in principle de- 
termined by the equation detM[q, w(q)] = 0. To make the 
result more physical, we decompose the complex fluctuation 
field ^{x) into its amplitude mode A(x) and phase mode 6{x), 
4>{x) - A(x) + i6{x). Then the fluctuation part of the efifective 
action takes the form 



c(2) 



i ^ ( r(Q) 6'(Q) )N(0 ( ^Jg j , (1 14) 



where the matrix N{Q) is defined as 
M+ + Mi2 



N(0 = 2 



-M: 



Ml 



11 

Mi2 



Here the quantities M*j are defined as 

Mf 1 (q, w) = i [Ml 1 (q, w) + Ml 1 (q, -w)] . 



(115) 



(116) 



We notice that M|j and Mj^j are even and odd functions of w, 
respectively. 

From the explicit form of Mi liQ), we have Mj^j (q, 0) - 0. 
Therefore the amplitude and phase modes decouple com- 
pletely at w = 0. Furthermore, using the saddle point condi- 
tion for the order parameter Aq, we find M+ (0, 0) = Mi2(0, 0), 
which ensures that the phase mode at q = is gapless, i.e., the 
Goldstone mode. 

We now determine the velocity Cj of the Goldstone mode, 
w(q) = Cjiql for w, |q| <s; minkjis^). For this purpose, we 
make a small q and oj expansion of N(2), 

M+ + Mi2 = A + C|qp - DcLT + ■■■ , 
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Here we note that the coefficient Q is proportional to the su- 
perfluid density and the superfluid phase stiffness J^- The 
explicit form of A, B, D, R and Q can be calculated as 
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The Goldstone mode velocity or the so-called sound veloc- 
ity in the superfluid state is given by 



Q 

B^/A + R' 



(119) 



The corresponding eigenvector of N is (A, 0) — (-ic\q\B/A, 1), 
which is a pure phase mode at q = but has an admixture 
of the amplitude mode controlled by B at finite q. Another 
massive mode, or the so-called Anderson-Higgs mode, has a 
mass gap 
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(A) Analytical Results for Large SOC. In the rashbon 
BEC limit A/kp » 1, we have ^ £6/2 and Aq <«; \iu\. There- 
fore, the coeflicients A, B, D, R and Q can be weU approxi- 
mated as 
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In this case, we find that B'/A » and therefore the ampli- 
tude and phase modes are strongly coupled. 

The sound velocity C j and the mass gap Mg read 



c, = 




(122) 



Using the relation d/a- - AjraQB/ms, the sound velocity re- 
covers the result for a weakly interacting rashbon gas, 



47raBBnB 



(123) 



Therefore, in the BEC limit, the quantity Cj/co depends only 
on the dimensionless parameter k - l/(Aas), where cq = 
ypjmjA. Using the results for otb and abb , we obtain 
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(124) 



The numerical result for the quantity c/cq is shown in Fig. 
|9] We find that it has a maximum near the point /c = 0, at 
K - -0.18. For the case Acu — » oo, we have 



Cs(Aas oo) - cq 



7 



\ 3(4 + V2) 



= 0.66co. 



(125) 
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FIG. 9: The sound velocity Cj of the Goldstone mode (divided by 
Co) in the RBEC regime as a function of the dimensionless parameter 
K= 1/(M). 



Using the expressions for a and d, we obtain the explicit 
form of the mass gap M„, 



4£b(£b 



— A 



J + 2 



(126) 



Therefore, in the BEC limit, the quantity Mg / depends only 
on the dimensionless parameter k - l/{Aas)- The numerical 
result is shown in Fig. [TO] We find that it is very small in the 
limit K -oo, and increases rapidly in the regime k > 0. For 
the case Aa^ — > oo, we have = 2A~ and therefore 
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FIG. 1 1 : The velocity of the Goldstone mode (sound velocity) Cj 
(divided by vp) as a function of A/kp for various values of 1/A;pa,. 
The red dashed lines corresponds to the analytical result (124). 



FIG. 10: The mass gap Mg of the Anderson-Higgs mode (divided by 
A^) in the RBEC regime as a function of the dimensionless parameter 
K = l/(Aas). 



(B) Numerical Results. Using the same trick in Section IV, 
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FIG. 12: The mass gap Mg of the Anderson-Higgs mode (divided 
by €p) as a function of A/kp for various va;lues of l/(fcpa,). The red 
dashed lines corresponds to the analytical result (126). 



we obtain 
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dimensionless quantities A,B,D,R and 2 are de- 
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Therefore, we have 



Q 



' B^-IA + R 



and 
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(131) 



where fp = kelm (m - 1) is the Fermi velocity for the non- 
interacting Fermi gas in the absence of SOC. 

Using the solutions of x\ and xz from the gap and num- 
ber equations, we can calculate the quantity Cj/yp and Mg/ep 
for given values of l/(A:pfl5) and Ajk^. The numerical results 
are shown in Fig. [TT|and Fig. [12] For large negative values 
of 1 / (kpOs) and A/kp — » 0, we recover the well-known result 
Cs - i^p/ V3 for weak coupling fermionic superfluids For 
negative values or small positive values of l/{kpas), the nu- 
merical result becomes already in good agreement with the 
analytical results (124) and (126) at A/kp ~ 1, which is con- 
sistent with the observation that the system enters the rashbon 
BEC regime at /l/A:p ~ 1. For large positive values of l/(^pfl 5), 
the numerical results are in good agreement with the analyti- 
cal results for all values of A/kp. 
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For very large Ajk-p, we find that the numerical results fit 
very well with the following scaling behavior 

±-.,mM±Y\ ^=4(if, (132, 

for both negative and positive values of 1 / (kpOs), as indicated 
from the analytical observations. 

VII. SUMMARY 

In summary, we have presented a comprehensive study of 
the BCS-BEC crossover problem in 3D Fermi gases with a 
spherical spin-orbit coupling which can be realized by a 3D 
symmetrical configuration of the synthetic SU(2) gauge field. 
The two-body problem, the superfluid ground-state properties, 
and the behaviors of collective excitations are studied. Ana- 



lytical results and interesting universal behaviors for various 
physical quantities at large SOC are obtained. We notice that 
there has been experimental proposal for the realization of 3D 
spherical spin-orbit coupling for cold fermionic atoms ifJsll . 
Therefore, it is interesting to test our theoretical predictions 
in future experiments of cold Fermi gases with 3D spherical 
spin-orbit coupling. 
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